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1 .  INTRODUCTION 


Radars,  operated  at  the  Kwajalein  Missile  Range  (KMR)  at  Kwajalein 
Atoll,  Marshall  Islands  on  the  island  of  Roi-Namur,  collect  weather  data 
on  hydrometeor  scatter  before  and  after  reentry  missions  of  importance 
to  military  programs.  Two  of  the  radars,  AI.COR  and  TRADliX,  part  of  the 
Kiernan  Reentry  Measurement  Site  (KREMS)  facility  operated  hy  the  Army 
Kwajalein  Missile  Range  Directorate,  BMD  Systems  Command  under  the  direc¬ 
tion  of  MIT  Lincoln  Laboratory,  are  used  for  the  collection  of  the  data 
on  clouds  and  precipitation.  These  investigations,  conducted  by  the  Air 
Force  Geophysics  Laboratory  for  diagnostic  measurements  and  studies,  arc 
in  support  of  U.S.  Air  Force  SAMSO/ABRES  programs  for  the  assessment  of 
the  severity  of  projected  and  actual  environmental  conditions.  The  ALCOR 
radar  operates  at  C-band  (5664  GHz)  and  has  been  used  for  weather  scan¬ 
ning  and  backscatter  data  collection.  The  12-meter  antenna  is  housed  in 
a  radome.  The  TRADEX  radar  has  a  26  m  reflector  (no  radome)  and  oper¬ 
ates  at  S-band  (2951  GHz)  for  the  same  purposes  on  this  type  of  data 
collection.  The  emphasis  in  this  report  is  on  the  documentation  and 
definition  of  the  parameters,  assumptions  and  constants  used  in  the 
backscatter  mode  when  the  radars  are  operated  and  data  are  digitally 
recorded  for  analysis  and  interpretation  in  terms  of  hydrometeor  scatter. 

The  use  of  the  KREMS  radars  for  weather  data  collection  has  been 
previously  documented  (Lewis,  1978)  with  emphasis  on  overall  weather 
measurement  capabilities  and  on  a  particular  software  program  called 
"MOIST"  also  developed  for  AFGL  for  use  at  sites  in  the  real  time  assess¬ 
ment  of  moisture  content  and  weather  severity  conditions  in  the  tropo¬ 
sphere.  Recently,  however,  the  radar  parameters  have  been  scrutinized 
again  with  regard  to  definitions  and  an  assessment  of  errors  and  bias 
terms  involved  in  calibrating  the  radars  and  processing  the  data.  This 
report  serves  as  a  documentation  of  the  updated  results  and  should  prove 
useful  to  those  involved  in  the  quantitative  analysis  of  the  backscatter 
data. 

As  a  result  of  this  study,  the  scattering  relationships.  Section  2, 
have  been  reviewed  with  a  special  emphasis  on  signal  processing  correc¬ 
tions  and  radar  calibration  uncertainties  in  Section  3.  The  signal  pro¬ 
cessing  effects  which  potentially  can  produce  bias  in  the  measurements 
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arc  incorporated  into  the  formulation  as  correction  terms.  The  radar 
calibration  uncertainties  must  be  considered  an  estimate  of  the  experi¬ 
mental  error  of  measurement  with  random  occurrences.  These,  sources  of 
measurement  error  are  summarized  in  Table  1  in  terms  of  RMS  (To)  uncer¬ 
tainty  in  decibels  upon  the  measurement  of  backscattert.J  radar  cross 
section  and  effective  radar  reflectivity  (directly  proportional) . 
Finally,  a  section  on  data  interpretation  is  included  (Section  4)  to 
assist  in  the  selection  of  hydrometer  scatter  due  to  rain  using  the  dual 
polarization  capability  of  the  KREMS  radars. 
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TABLE  1 


SUMMARY  OF  SOURCES  OF  MEASUREMENT  UNCERTAINTY 
IN  KREMS  RADAR  HYDROMETEOR  BACKSCATTER  DATA  US I NO 
THE  RADAR  EQUATION  AND  IN  RADAR  CALIBRATION1, 


2 

Estimated  Maximum 

Uncertainty  (dB) 

Source  and  Affected  Term  * 

ALCOR 

TRADI'.X 

Quantization,  Pout 

0.25 

0.5 

Attenuators,  A 

0.05 

0.05 

4 

Range,  R 

0.1 

0.1 

Transmitted  Power,  P 

0.1 

0.1 

Beamwidths,  064 

0.5 

0.5 

Calibration  Constant,  KRCS 

0.6 

0.6 

Receiver/Matched  Filter 
Stability  (short  term),  PR 

0.5 

0.5 

Propagation  Effects  and 

System  Loss  Differences 
(calibration-to-test)^* ,  L 

0.5 

0.3 

Total  Estimated  Maximum 
Uncertainty 

±2.6  dB 

+2.7  dB 

Estimated  RMS  (lo) 

Uncertainty 

±1.0  dB 

±1.0  dB 

1,see  section  3.2 
2 

*see  equations  (15)  and  (16)  and  associated  relationships 


'"omitted  from  treatment  are  possible  polarization  mis-match  losses,  m, 
and  radome  attenuation  changes  both  assumed  here  to  have  zero  dB 
uncertainty  throughout  a  hydromoteor  measurement 


J 
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2.  EFFECTIVE  REFLECTIVITY  FOR 
DISTRIBUTED  HYDROMETEOR  SCATTER 


2.1  Backscatter  from  Rain 


For  backscatter  from  rain  drops,  the  radar  reflectivity  or  scatter¬ 
ing  cross  section  per  unit  volume,  8,  is  the  volume  summation  of  the 
reflectivities,  8.,  from  individual  scatterers: 


B  = 


1_ 

V 


total  radar  scattering  cross  section 
volume  of  scattering 


For  radar  wavelengths  large  compared  to  the  drop  size  such  that  Rayleigh 
scattering  holds,  i.e.  for  small  spherical  water  droplets  and  radar  fre¬ 
quencies  below  7  GHz;  the  particle  backscatter  cross  section  is: 

°i  3  4  |i(|2  Di6  »  lKl  3  (1) 

where  c  is  the  dielectric  constant  for  the  drop; 

A  is  the  radar  wavelength;  and 
D.  is  the  diameter  of  the  drop. 


For  volumetric  scatter,  the  total  cross  section  is: 


I 


|K|2  I  D,6 

voi 


O) 


2 


(2) 


assuming  no  polarization  mismatch  in  electromagnetic  scattering. 

fn  general,  the  scattering  cross  section  per  unit  volume  is  expressed 
(Grnne,  1973;  also  1974)  for  radar  scatter  as: 

r  =  I_  |K|2  m  Z  (m2/m3)  (3) 

A 

where  m  =  polarization  mismatch  factor;  and 

Z  =  radar  reflectivity  factor  proportional  to  the 

sixth  power  of  particle  diameters  {Z  =  -y  £  D.^}. 

vof.  1 

Usually,  the  radar  reflectivity  is  expressed  as  the  effective  radar 
reflectivity  for  rain,  2^,  conventionally  in  units  of  matin''  where  Zg  is 
defined  as: 
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with  | K( 2  =  0.93,  (e  n.  81  for  water). 


The  elemental  volume,  V,  illuminated  by  a  uniform  beam  and  a  rec¬ 
tangular  pulse  in  a  single  radar  resolution  cell  is: 

V  =  ir(R  |  )  (R  f  )  c  J  (5) 

where  R  =  slant  range  to  scattering  volume; 

0,<j>  =  radar  angular  beamwidths; 
t  =  pulse  width;  and 
c  =  velocity  of  light. 

The  radar  equation  for  backscatterc J  received  power  in  a  resolution 
cell,  assuming  uniform  B  over  the  volumetric  element  is: 


where  P  =  transmitted  power  delivered  to  antenna; 

G  =  antenna  gain;  and 

I,  =  total  system  and  propagation  losses  (two-way) , 

(I.  >  1). 

A  more  accurate  representation  of  the  actual  scattering  volume  is 
obtained  by  approximating  the  antenna  pattern  as  Gaussian  shaped  in  the 
elevation  and  orthogonal  angle  planes  (Probert -.Jones,  1962),  and  the 
time-domain  response  as  a  Gaussian  shaped  pulse.  The  resulting  repre¬ 
sentation  for  the  volume  element  of  scattering  with  a  single  resolution 
cell  is  given  by: 


where  0,$  =  'j  power  beamwidths  of  antenna  (one-way) ; 

ei/2  =  I)q-.  /w/4  In  2  =  1.064  1)^,  (see  eqn.  (13)); 
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I>0  =  'i  power  range  resolution  »  ctq/2;  and 
To  =  1,2  Power  Pulsc  width  (compressed). 


f4). 


The 

(6) 


average  power  per  resolution  cell  may  be  expressed  by  combining 
and  (7)  as: 


P 


r 


v>2  y.s 

(4ir)  Vl, 


C«y 


The  scattered  power  in  terms  of  radar  cross  section,  cr,  from  an 
equivalent  point  scatterer  at  the  same  position  as  the  elemental  volume 


where  a  -  radar  cross  section  (m  ). 


liquating  (8)  and  (9)  permits  the  equivalent  radar  cross  section  of 
the  scattering  cell,  a,  to  be  expressed  in  terms  of  effective  radar 
reflectivity,  2^: 


<» 


,6|k|2„  ze-io-‘Ve«  i>„7^ 

\ 4  8  ln2 


(10) 


2.2  Backscatter  Reflectivity  for  General  Hydrometeors 


Given  that  the  radar  cross  section,  cr,  is  a  measurable  quantity 
derived  from  the  calibrated  received  power  output  of  the  radar  (the 
other  terms  and  radar  constants  of  (9)  arc  measured  or  known),  the 
effective  radar  reflectivity  from  a  scattering  cell  may  be  obtained  by 
solving  equation  (10)  for  Z^: 


1018  A4  8  ln2 


ir  y/’mlK'  |  ~o<j> 


(mm^/m3) 


(ID 


4  ln2 


where  general  hydromctcor  scatter  is  relatable  to  the  equivalent  for 
water  droplets: 


|K'|  =  | c- 1 /r+2 |  •  1/p; 


c  =  dielectric  constant  for  the  type  of  hydrometeor 
producing  effective  scatter;  and 
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p  a  ratio  of  density  of  the  hydromcteor  to 
the  density  of  water. 


Equation  (11)  may  alternatively  be  expressed  as: 


a  C 

7  * 


or  dBZ  =*  10  log  Zg  =  a(dBSM)  -  20  log  R  +  10  log(C/F) 


where 


C  =  calibration  constant 


{r  is  expressed  in  meters}; 


a  *  radar  cross  section  expressed  in  decibels 
relative  to  a  square  meter  (dBSM) ;  and 

F  =  statistical  processing  correction  factor 

which  depends  on  signal  processing  algorithms 
in  use  with  radar  data.  _ ... 


When  backscatter  is  from  hydrometeors  other  than  rain  from  this 
formulation,  it  is  necessary  to  use  the  appropriate  value  of  K*  for  the 
type  of  scatter  involved.  In  the  MOIST  program  processing  the  height  of 
the  scattering  volume  is  obtained  from  refraction  corrected  radar  para¬ 
meters  in  real  time  and  used  as  a  criterion  on  which  to  base  the  use  of 
a  K'  for  ice  versus  a  K'  for  water.  The  height  used  in  the  program 
(selectable)  is  that  of  the  melting  layer  (^  0®C  isotherm  height)  and  is 
typically  4.6  km  for  Kwajalein.  This  altitude  corresponds  to  the  lower 
edge  of  the  radar  "bright-band"  below  which  melting  takes  place. 
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3.  CALIBRATION  OF  ALCOR  AND  TRADEX  RADARS 


The  ALCOR  (C-band)  and  TRADEX  (S-band)  radars  operate  in  a  narrow 
band  pulse  compression  mode  (CHIRP)  for  hydrometeor  scatter  measurements. 
The  compressed  output  pulse  (matched  filter  response)  may  be  modeled  as 
approximately  Gaussian  shaped.  This  assumption  is  consistent  with  quoted 
3-  and  6-dB  pulse  widths  and  the  shape  of  time  response  plots  obtained 
during  calibration  runs.  For  the  radars,  the  sampling  volume  is  defined 
by  the  antenna  beamwidth  and  the  response  of  the  pulse  compression  cir¬ 
cuits.  In  angle,  the  volume  definition  is  given  by  Probert-Jones  (1962) 
as  indicated  in  the  previous  section.  In  range,  the  volume  is  defined 
by  the  integral  over  the  time  response: 


cr 

2 


-4  In2(x/D  ) 
e 


d* 3  DoVrfer 


1.064  D 

o 


(13) 


CT 

where  Dq  =*  — —  =  pulse  resolution  width  (tq  =>  3  dB  pulse  width) . 

From  the  derived  expressions  (11)  and  (12)  the  effective  radar  reflec- 
6  3 

tivity  (mm  /m  )  may  be  obtained  from  the  measured  radar  scattering  cross 
section,  range  to  the  scattering  volume,  calibration  constant  C,  and 
processing  factor  F.  The  radar  parameters  and  calibration  constants  are 
given  in  Table  2  for  the  current  best  estimate  values  for  the  two  KREMS 
radars. 


3.1  Signal  Processing  Corrections 

The  processing  correction  factor,  F,  is  used  to  adjust  for  the 
biases  introduced  in  the  average  value  of  the  randomly  varying  signal 
received  from  an  ensemble  of  point  scatterers  within  the  scattering 
volume.  The  processing  correction  factor  depends  upon  the  particular 
digital  signal  processing  scheme  employed  in  the  data  reduction  chain  to 
obtain  the  equivalent  scattered  cross  section  used  in  equation  (12). 

The  correction  factor  accounts  for  biases  due  to:  digital  quantization 
(sometimes  referred  to  as  truncation  error);  detection  cell  selection 
algorithm  biases;  and  biases  associated  with  averaging  the  logarithm 
of  the  amplitude  rather  than  the  value  of  the  received  power.  In  contrast 


TABLE  2 


RADAR  PARAMETERS*  FOR  VOLUMETRIC 
SCATTER  FROM  HYDROMETEORS 


Parameter 

ALCOR 

TRADEX 

Units 

A 

0.05292 

0.1016 

(m) 

m 

1.0 

1.0 

- 

0  =  4 

5.3  x  10-3 

5.0  x  10‘3 

(rad.) 

D 

o 

37.5  m 

10.4 

(m) 

|K’|2  ice 

0.209 

0.209 

- 

|K’ | 2  water 

0.933 

0.934 

- 

C  ice 

82.9  (142.9) 

100.3  (160.3) 

(dB) 

C  water 

76.4  (136.4) 

93.8  (153.8) 

(dB) 

Values  in  parentheses  are  for  range  expressed  in  meters, 
R  is  in  km 

otherwise 

to  most  weather  radars,  the  KRF.MS  radars  are  calibrated  effectively  using 
the  center  of  the  quantization  unit.  A  quantization  correction,  there¬ 
fore,  is  not  required  since  no  bias  should  be  introduced.  The  detection 
cell  selection  algorithms  and  received  pulse  averaging  schemes  introduce 
significant  amplitude  biases  which  may  be  corrected. 

The  two  statistical  sources  of  bias  produced  by  the  processing 
scheme  used  with  the  KREMS  radars  are  those  introduced  by  averaging  the 
logarithm  of  the  received  signal  (both  radars)  and  the  "greatest-of-four" 
(consecutive  range  samples)  algorithm  used  in  the  real-time  MOIST**  pro¬ 
cessing  of  the  ALCOR  data.  The  probability  density  function  of  the 
average  of  a  number  of  independent  samples  of  received  power  in  a  reso¬ 
lution  cell  approaches  a  Gaussian  function  for  a  large  sample  size 
(exceeding  10;  see  Smith,  1964).  For  the  case  in  which  the  receiver  has 
a  logarithmic  amplitude  response  (as  for  the  KREMS  radars),  the  average 

*Note  that  the  parameters  in  Table  2  have  been  revised  from  those  pub¬ 
lished  by  Lewis  (1978),  mostly  due  to  improved  radar  beamwidth  and  range 
resolution  information. 

**M0IST  is  a  software  program  used  at  KREMS  for  real-time  estimation  of 
liquid  water  content  using  either  ALCOR  or  TRADEX  radar  data  output. 
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of  the  log  power  (10  log  PR)  is  obtained.  The  desired  quantity  10  log  PR 
(logarithm  of  the  average  power)  is  obtainable  with  a  correction  term: 

10  log  PR  =  10  log  PR  +  2. SI  dB  (14) 

The  bias  correction  term  in  F  due  to  averaging  log  power  is  therefore 
2. SI  dB. 

For  ALCOR,  with  the  data  obtained  in  real-time  and  processed  through 
the  MOIST  program,  the  greatest-of-four  consecutive  range  samples  are 
selected  prior  to  averaging  for  post  detection  processing.  The  samples 
are  not  independent  since  the  half-power  range  resolution  is  37.5  meters 
but  samples  are  taken  every  IS  meters.  The  correction  for  the  effect  of 
processing  using  the  greatest-of-four  detection  algorithm  can  be  esti¬ 
mated  either  through  simulation  (Smith,  1964)  or  by  direct  measurement. 
The  results  in  Table  3  were  obtained  by  direct  measurement  by  employing 
the  MOIST  processing  scheme  off-line  using  data  obtained  from  an  ALCOR 
data  tape.  A  summary  of  the  combined  processing  factor  F  for  data  pro¬ 
cessed  with  the  radars  is  included  in  Table  3  together  with  the  corrected 
calibration  constant,  C/F. 

3.2  Radar  Cross  Section  Calibration 

The  radar  parameters  contained  in  parentheses  in  equation  (9)  must 
either  be  measured  or  deduced  to  obtain  a  correct  estimate  for  radar 
cross  section,  a,  given  the  measured  value  of  received  power.  The  KREMS 
radars  are  calibrated  by  tracking  a  calibration  sphere  of  known  radar 
cross  section.  The  constants  with  the  exception  of  transmitted  power  and 
range  are  deduced  from  the  radar  tracks  of  several  spheres.  The  result¬ 
ing  expression  for  an  unknown  scattering  cross  section  is: 


where  R  =  measured  range  to  scatterer  (meters) ; 

^RCS  =  const;mt  obtained  from  external  sphere  target 
RCS  calibration 

=  ^  "90  dB  for  ALCOR),  (m-2); 

GT 

PR  =  measured  received  power  (watts) 
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TABLE  3 


+  A 


(dB); 


3  P 

OUT 

PqUT  3  power  received  after  attenuation  (watts); 

A  3  total  attenuation  in  receiving  system  (dB); 

L  3  total  two-way  system  loss  including  transmission 
lines  and  propagation  (L  >  1); 

G  »  on-axis  radar  antenna  transmit  and  receive  gain 
factor;  and 

A  a  radar  operating  wavelength  (meters) . 

Expressed  in  decibels,  equation  (15)  becomes: 

a(dBSM)  =  poirr  +  A  '  PT(dBW)  +  40  log  R  *  kkcs 
where  ^OUT  =  Table  A/D  reference  (dBW)  +  Table  A/D  difference  (dB) 

Figure  1  shows  a  simplified  block  diagram  of  the  ALCOR  radar,  for 
example  (narrow  band  waveform  configuration) ,  used  in  most  weather  radar 
applications.  The  receiver  chain  is  calibrated  by  signal  injection  pre¬ 
ceding  the  logarithmic  detector.  The  output  is  A/D  converted  prior  to 
digital  recording.  The  characteristic  curve,  A/D  output  versus  input 
level,  may  be  obtained  by  either  a  curve-fitting  or  a  point-by-point 
calibration  at  each  quantization  level.  The  actual  procedure  with  ALCOR 
is  to  produce  a  calibration  table  (dB  referenced  to  the  highest  A/D 
level,  127)  for  each  A/D  count  (level)  over  the  response  dynamic  range. 

The  output  power  for  equation  (16)  is  therefore  obtained  as  the  table 
A/D  difference  value  (dB)  added  to  the  table  A/D  reference  value  (dBW) 
obtained  for  the  reference  count  level.  This  look-up  table  method 
results  in  a  calibration  characteristic  for  the  log  detector  and  A/D 
chain  which  has  no  net  quantization  bias  but  varies  about  the  mean  fitted 
curve  by  a  residual  amount  typically  less  than  one  quantization  unit  in 
magnitude  (0.55  dB  for  ALCOR,  1.11  dB  for  TRADEX) .  There  is,  however, 
a  maximum  uncertainty  associated  with  this  procedure  of  the  order  of  'i 
of  the  quant i zat ion  value  (%  0.25  and  %  0.5  dB,  respectively),  depending 
on  the  instantaneous  signal  level  within  the  dynamic  range.  Assuming 
some  signal  fluctuation  throughout  the  period  of  averaging,  the  uncer¬ 
tainty  associated  with  this  measurement  procedure  will  be  reduced  below 
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Figure  1  Simplified  block  diagram  of  the  ALCOR  radar  used  for  weather  radar  applications 


that  expected  on  a  single  pulse  basis.  The  uncertainty  (one  standard 
deviation)  is  equal  to  the  magnitude  of  the  quantization  step  divided  by 
the  square  root  of  12  times  the  number  of  independent  samples.. 

The  other  terms  remaining  in  equation  (16)  may  also  have  uncertain¬ 
ties  associated  with  their  measurement.  Some  of  these  are  negligible 
such  as  those  produced  by  the  attenuator  calibration  (less  than  +0.05  dB 
error)  and  that  of  range  R  estimation  (less  than  +0.1  dB)..  The  trans¬ 
mitted  power,  P^,  is  monitored  and  factored  into  the  calibration  proce¬ 
dure  by  the  radar  real-time  program.  Though  the  transmitted  power  level 
may  vary  slowly,  or  on  a  miss ion- by-mission  basis,  the  calibration  proce¬ 
dure  removes  all  but  a  random  component  which  is  estimated  to  be  less 
than  +0.1  dB.  The  term  which  has  the  largest  uncertainty  is  the  calibra¬ 
tion  constant  usually  obtained  from  sphere  calibrations  preceding 

or  following  a  mission  run.  The  uncertainty  associated  with:  a  sphere  of 
’’known'*  radar  cross  section  (RCS)  depends  on  the  size  of  the  sphere  rela¬ 
tive  to  the  radar  wavelength.  For  scattering  from  a  sphere  within  the 
optic  regime  (ka  >  10,  k  2-ir/X ,  a  =  sphere  radius)  the  RCS  uncertainty 
is  less  than  +0.5  dB  assuming  a  perfectly  spherical  shape.  The  20- inch 
diameter  spheres,  flown  for  radar  calibration  before  and  after  missions 
at  Kwajalein,  are  clearly  within  the  optical  criterion  cited,  at  the  radar 
wavelengths  (ka  -v  30  for  ALCOR  C-band,  ka  %  16  for  TRADEX  S-band) . 

Other  factors  may  enter  into  the  procedure  of  calibration  beyond  the 
uncertainty  of  knowing  the  precise  sphere  RCS.  Since  K  contains 
propagation  losses  as  well  as  system  parameters  there  may  be  differences 
between  those  encountered  during  sphere  calibration  and  those  during 
actual  operation  for  hydrometeor  scatter.  Some  of  these  are:  different 
path  attenuations  (due  to  absorption  or  hydrometeor  attenuation)  associ¬ 
ated  with  the  different  atmospheric  paths;  differences  associated  with 
tropospheric  refraction  effects;  variations  in  radome  attenuation  due  to 
rain  at  AI.COR;  uncertainties  associated  with  target  positioning  within 
the  radar  beams;  and  gradual  changes  in  system  gains,  losses  or  matched 
filter  response  from  the  time  of  calibration  to  the  time  of  measurements. 
Generally,  the  differences  are  randomly  occurring  events,  the  magnitude 
of  which  can  only  be  estimated  through  many  observations.  Typically, 
the  changes  in  calibration  constant,  an  indication  of  some  of  these 

effects  occurring,  are  not  more  titan  +l-to-2  dB  from  one  sphere  calibration 


14 


/ 

/ 


/ 


to  the  next  and  this  principally  charged  to  variation  in  channel  gain. 
During  an  operational  period,  the  radar  system  itself  is  considered 
stable  to  about  ±0.5  dB.  The  uncertainty,  therefore,  associated  with 
depends  on  the  proximity  in  time  and  space  of  the  scatter  measure¬ 
ment  to  the  sphere  calibration  and  upon  operator  care  in  system  set-up. 

It  is  estimated  that  errors  in  determining  KRCS  are  ±0.6  dB  (±1.8  dB  at 
3o) including  the  sphere  RCS  uncertainty  when  conditions  are  nominal,  i.e. 
no  rain  occurrs  at  the  radars  and  sphere  calibrations  are  made  within  a 
few  hours  of  the  measurement  run.  This  estimate  is  based  upon  an  analy¬ 
sis  of  sphere  calibration  results  for  TRADEX  during  a  four-month  period 
in  1977  which  resulted  in  an  RMS  error  in  sphere  calibrations  to  be 
within  1  dB  (see  also  Crane,  1978).  The  errors  for  a  particular  mission 
are  less  than  the  variability  for  a  number  of  widely  spaced  sphere  cali¬ 
brations.  Errors  associated  with  other  propagation  effect  differences, 
e.g. ,  gaseous  attenuation,  hydrometeor  attenuation  and  system  loss  changes 
(from  calibration  to  mission),  are  estimated  at  0.5  dB  for  ALCOR  and 
0.3  dB  for  TRADEX.  The  uncertainties  are  summarized  in  Table  1. 


4.  HYDROMETEOR  SCATTER  DATA  AND  INTERPRETATION 


Radar  measurements  of  hydrometeor  scatter  at  Kwajalein  are  collected 
for  diagnostic  as  well  as  operational  purposes.  When  diagnostic  measure¬ 
ments  are  taken,  an  instrumented  aircraft  is  flown  to  sample  the  precipi¬ 
tation  and  make  direct  measurements.  The  radars  are  slaved  to  continu¬ 
ously  track  a  scattering  volume  ahead  of  the  aircraft  path  such  that  data 
may  be  correlated  between  the  two  sets  of  measurements.  An  arbitrary 
sample  of  data  from  a  correlation  run  is  examined  next,  simply  to  illus¬ 
trate  the  type  of  data  collected  and  to  provide  a  basis  for  discussing 
the  interpretation  of  rain  scatter  data  in  particular. 

4.1  Short  Sample  of  Rain  Scatter  Data 

An  example  of  rain  scatter  measurements  is  shown  in  Figure  2  for  one 
of  several  correlation  runs  (taken  near  Roi-Namur  Island)  on  June  26,  1979. 
This  short  selected  sample  occurred  for  a  data  collection  pass  through  a 
light  rein  shower  region  which  was  relatively  inhomogeneous.  The  reflec¬ 
tivity  data  are  from  MOIST  data  processing  at  a  rate  of  one  averaged 
value  per  second.  Both  ALCOR  and  TRADEX  radar  data  are  shown  on  the 
figure  for  comparison.  Typically,  the  effective  radar  reflectivity  in 
decibels,  dBZ  (varying  in  this  case  from  -20  to  +30  dBZ),  correlates 
well  between  measurements  taken  by  the  two  radars  with  similar  scatter¬ 
ing  volume  resolution  (as  noted  in  Table  2).  Shown  also  in  the  figure  is 
the  ratio  of  principal  (PP)  to  opposite  (OP)  received  backscatter  polar¬ 
ization.  The  principal  polarization  for  a  circularly  polarized  radar  is 
defined  as  left-circular  when  right-circular  polarized  waves  are  trans¬ 
mitted  (opposite,  is  right-circular). 

4.2  Use  of  Cross  Polarization  as  an  Editing  Criterion 

When  operating  under  conditions  where  rain  is  the  supposed  cause  of 
hydrometeor  scatter,  one  can  use  the  radar  cross  polarized  channel  as  an 
indication  or  gage  of  the  nature  of  the  radar  backscatter  as  rain  rather 
than  another  form  of  hydrometeor.  The  scattered  components  of  the 
electromagnetic  field  received  by  the  radar  antenna  as  cross  polarized, 
will  increase  (in  the  case  of  backscatter)  as  the  scatterers  depart  from 
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Example  of  effective  radar  reflectivity  measurement,  dBZ, 
obtained  for  rain  on  AI.COR  and  TRADEX  radars  (variable 
conditions  near  island  of  koi-Namur) 


having  a  spherical  shape.  Depolarization,  or  a  rise  in  the  cross  polar¬ 
ized  channel  signal  power,  is  an  indication  of  the  non-symmetrical  nature 
of  the  scatterers  such  as  would  be  expected  from  prolate  or  oblate  spher¬ 
oids  when  viewed  at  an  arbitrary  angle.  Spherical  rain  drops,  in  con¬ 
trast,  should  backscatter  no  field  component  with  the  polarity  opposite 
to  that  defined  as  principal  on  the  axis  of  the  radar  beam.  The  ratio  of 
principal  to  orthogonal  PP/OP  polarized  returns  does  not  go  to  infinity 
however  due  to  channel  isolation  limits  of  the  antenna  system  (approxi¬ 
mately  26  to  28  dB) . 

The  oblateness  of  raindrops  is  a  function  of  size  and  hence  rain 
rate.  Since  these  characteristics  arc  measurable  and  known  as  a  function 
of  rain  temperature,  one  can  derive  a  relationship  for  the  polarization 
ratio  terms  of  rain  rate  and  total  scattered  reflectivity.  The  ratio  of 
power  received  in  the  principal  to  the  power  received  in  the  opposite 
channel  for  non-canted  raindrops  at  a  radar  frequency,  F,  is  therefore: 


PP/OP  = 


P-1 

P+1 


2 


(17) 


where  P  -  [Zy/Z^S 

Zy  =  reflectivity  for  vertical  linear  transmitted  and 
received  polarization  at  frequency  F;  and 

Zjj  =  reflectivity  for  horizontal  linear  transmitted 
and  received  polarization  at  frequency  F. 


F.xprcssed  in  decibels,  (17)  is: 


PP  -  OP  (dB)  =  20  log1Q 


P-1 

P+1 


(18) 


Once  again,  this  result  is  onl>  crue  in  the  ideal  sense  for  oblate 
spheroidal  raindrops,  with  a  vertical  orientation  of  the  axis  of  symmetry, 
scattering  on-axis  to  the  radar  beam,  with  a  zero  depth  of  penetration 
into  the  scattering  volume,  and  received  on  a  perfectly  polarization 
balanced  radar  antenna  system  at  very  high  signal  to  noise  ratios.  With 
a  non-ideal  real  radar  system,  some  allowance  must  he  made  for  degradation 
due  to  each  of  these  sources  of  departure.  It  is  still  of  utility,  how¬ 
ever,  to  derive  a  rule  of  thumb  for  judging  the  quality  of  radar  rain 
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backscatter  by  establishing  an  editing  criterion  albeit  somewhat  arbi¬ 
trary  (3  dB  or  more  departure  from  ideal,  as  a  function  of  dBZ  level), 
without  a  careful  investigation  into  each  potential  source  of  polariza¬ 
tion  degradation.  As  a  result,  the  included  table  of  minimum  editing 
rules.  Table  4,  may  be  considered  usable  as  an  estimate  of  the  necessary, 
but  not  sufficient,  conditions  to  select  the  total  backscattcrcd  dBZ 
levels  as  truly  rain  scatter  in  nature.  The  table  is  for  rain  at  1 0° C 
and  at  a  radar  frequency  of  4  GHz  (approximately  mid  range  between  the 
ALCOR  and  TRADEX  radar  frequencies).  Quite  obviously,  further  study 
into  this  type  of  classification  could  perhaps  yield  a  refinement  of  the 
criterion  and  a  similar  means  of  classification  f's  other  habits  of 
hydrometeor  activity  as  well  as  for  rain. 

TABLE  4 

EDITING  RULES  FOR  ASSISTING  IN  THE 
SELECTION  OF  USABLI-  RADAR 
BACKSCATTER  DATA  FROM  RAIN 

Range  of  dBZ  Minimum  Rule*  Ideal  (Theory) 


If  dBZ  <  0 

9 

then 

PP-0P 

> 

20 

dB 

off  scale 

0  £  dBZ  < 

10, 

then 

PP-0P 

> 

19 

dB 

- 

10  <  dBZ  < 

20, 

then 

PP-OP 

> 

18 

dB 

30 

20  <  dBZ  < 

30, 

then 

PP-OP 

> 

17 

dB 

25 

30  <_  dBZ  < 

40, 

then 

PP-OP 

> 

15 

dB 

20 

40  £  dBZ  < 

50, 

then 

PP-OP 

> 

13 

dB 

17 

50  <  dBZ  < 

60, 

then 

PP-OP 

> 

12 

dB 

15 

60  <  dBZ, 

then 

PP-OP 

> 

11 

dB 

14 

’Somewhat  arbitrary  at  lower  dBZ  values  due  to  lack  of  knowledge  of 
noise  statistics  influencing  the  PP/OP  ratio 


5.  SUMMARY  AND  CONCLUSIONS 


The  ALCOR  (C-band)  and  TRADEX  (S-band)  radars  at  Kwajalein  are  well- 
calibrated  instrumentation  radars  which  serve  among  many  other  roles  to 
provide  measurements  of  the  environmental  conditions  of  hydrometeor 
activity  within  the  reentry  corridor  of  the  Kwajalein  Missile  Range  (KMR). 
This  report  defines  the  essential  radar  parameters,  assesses  significant 
sources  of  error  in  these  parameters  and  documents  the  most  recent  best 
estimate  of  the  appropriate  radar  constants  to  vse  in  data  processing. 

A  discussion  of  the  calibration  procedures  is  also  given  with  a  recogni¬ 
tion  of  the  ultimate  importance  of  sphere  calibrations  in  close  space  and 
time  proximity  to  the  hydrometeor  scatter  measurements  of  interest  for 
minimizing  the  calibration  uncertainties.  The  overall  RMS  (la)  uncer¬ 
tainties  associated  with  the  radar  measurement  are  estimated  at  1.0  dB 
each  for  ALCOR  and  TRADEX  as  summarized  in  Table  1.  This  assumes  that 
the  bias  terms  associated  with  signal  processing  are  appropriately 
accounted  for  by  the  use  of  the  best  estimate  correction  factor  given 
in  Table  3  and  that  the  radar  parameters  given  in  Table  2  do  not  change. 
Implicit  in  this  analysis  is  the  assumption  that  a  zero  dB  polarization 
mis-match  loss  occurs  for  hydrometeor  scatter  and  that  excess  losses, 
such  as  radome  attenuation  changes  due  to  rain  or  condensation  of  mois¬ 
ture  on  the  ALCOR  radome,  do  not  occur  during  the  measurements. 

An  initial  attempt  to  classify  rain  scatter  on  the  basis  of  prin¬ 
cipal  to  opposite  polarization  ratio  is  presented.  This  tool,  useful 
for  data  editing  and  interpretation,  suggests  that  further  means  of 
classifying  the  hydrometeor  backscatter  may  be  intrinsically  available 
in  the  different  radar  polarizations  and  worthwhile  for  further  investi¬ 
gation  at  Kwajalein. 
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